GRAIN ALIGNMENT IN THE TAURUS DARK CLOUD 
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ABSTRACT 



Variations in polarization efficiency (p/A) of interstellar grains as a function of environment 
place vital constraints on models for the mechanism of their alignment. We use polarimetry 
of background field stars to investigate alignment in the Taurus dark cloud as a function 
of optical depth for extinctions of Ax < 2.5 (Ay < 25). Results show a strong systematic 
trend in polarization efficiency with extinction. We discuss the possibility of a relationship 
between p/A and icy grain mantle growth, since mantles are detected in Taurus only for 
extinctions where p/A is lowest. A simple density model of alignment by paramagnetic 
alignment (Davis & Greenstein 1951) predicts that the magnetic field strength (B) scales 
with density to the power n = 0.28. 

Subject headings: dust, extinction — polarization — ISM: magnetic fields — ISM: indi- 
vidual (Taurus Cloud) 
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1 INTRODUCTION 



The light from stars obscured by interstellar dust is almost always partially plane 
polarized at visible wavelengths, a phenomenon attributed to extinction by aspherical dust 
grains which are aligned to some degree by the interstellar magnetic field (see Whittet 
1995 for a review). It is widely agreed that alignment is the result of an interaction 
between the rotational dynamics of the grains and the ambient magnetic field, as originally 
proposed by Davis & Greenstein (1951). Paramagnetic relaxation results in the grains 
tending to become oriented with their angular momenta parallel (and hence their long axes 
perpendicular) to the magnetic field lines. Although this mechanism is qualitatively highly 
successful, a detailed quantitative theory of magnetic alignment that is consistent with all 
existing observational constraints is yet to be formulated. The efficiency of alignment 
depends on both the properties of the dust grains themselves and on the environment in 
which they exist (e.g. Jones & Spitzer 1967; Purcell 1979; Mathis 1986; DeGraff et al. 1995). 
Polarization studies, in conjunction with measurements of other physical and chemical 
characteristics of dense clouds, may thus be used as observational tools for constraining 
the alignment mechanism. 

The observed degree of linear polarization (p) correlates positively with the magni- 
tude of extinction (A) due to dust in the line of sight, but with scatter much greater 
than can be accounted for by observational errors alone (e.g. Serkowski et al. 1975). The 
polarization efficiency of the interstellar medium (ISM), represented by the ratio p/A at 
some wavelength, is highly non-uniform. In general, p/A falls below some upper bound- 
ary 3 representing optimum efficiency. Scatter in p/A below this optimum value might 
arise due to a number of factors. In a uniform medium with a uniform magnetic field, 
p/A varies systematically with the angle between the line of sight and the field direction, 
and is maximum when this angle is 90°. If the magnetic field is inhomogeneous, the di- 
rection of alignment may change along a given line of sight, leading to a reduction in net 
polarization. Similar depolarization effects may be observed if two or more clouds with dif- 
ferent field orientations exist within the same line of sight, and if the spectral dependence 
of polarization also changes from cloud to cloud, the degree of depolarization becomes a 
function of wavelength. Of greater interest are variations in p/A induced by changes in 
the efficiency of alignment as a function of physical conditions. In this paper, we attempt 
to isolate the effect of physical environment on polarization efficiency in order to improve 
our understanding of the alignment mechanism. 

The effects of viewing geometry and depolarization may be minimized by appropriate 
choice of sample. Our approach is to select data for stars obscured by an individual cloud 
displaying a relatively uniform magnetic field structure. The Taurus dark cloud is an 
excellent candidate. As it is relatively nearby (~ 140 pc, Elias 1978) and ~ 20° from the 
Galactic plane, the extinctions (and polarizations) of stars in this area of the sky arise 

3 In the diffuse ISM, a limiting value given by Pm&x/Ay < 3.0 % mag -1 is generally 
assumed (e.g. Serkowski et al. 1975), where p max is the peak percentage polarization in the 
Serkowski formulation, Ay is the visual extinction, and Ay = 3Eb-v has been assumed. 
Whittet et al. (1994) recently found values of Pm&x/Ay up to ~ 4.5 % mag -1 in the 
Chamaeleon I dark cloud. 
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almost entirely within the cloud itself (Straizys & Meistas 1980; Kenyon et al. 1994). The 
position angle of polarization is independent of wavelength in 15 out of 21 stars studied by 
Whittet et al. (1992), consistent with a general absence of discrete cloud components with 
different wavelength dependencies of polarization. The distribution of polarization vectors 
in the plane of the sky is fairly uniform over the entire cloud (Moneti et al. 1984; Tamura 
et al. 1987; Goodman et al. 1992). Although small-scale regions of magnetic complexity 
associated with individual condensations undoubtedly occur, the macroscopic field (as 
sampled by field stars) does appear to be generally lacking in major inhomogeneities. We 
therefore surmise that the degree of grain alignment is the primary factor determining the 
polarization efficiency, p/A } in the Taurus dark cloud. 

A further advantage of this choice is the fact that the Taurus cloud has been well 
studied in the previous literature and consequently its physical and chemical properties 
are reasonably well constrained. Studies of the Taurus- Auriga cloud complex show that 
it possesses a filamentary structure (e.g. Kleiner & Dickman 1984, Cernicharo et al. 1985, 
Wouterloot & Habing 1985), with small dense cores observed in line emission of molecules 
such as NH3 (e.g. Myers & Benson 1983) and CS (Zhou et al. 1989). A complex chemistry 
is evident in these cores, with such varied molecules as C3N (Friberg et al. 1980), HNCO 
(Brown 1981), CH3CCH (Irvine et al. 1981), and aromatic hydrocarbons (Hanner et al. 
1994) detected in the gas phase. They have typical masses of a few Mq, gas temperatures 
of ~ 10 K, and densities ~ 10 4 cm -3 (e.g. Wilson & Minn 1977), and are thought to be 
active sites of star formation (e.g. Lada et al. 1993), associated with T Tauri stars (Myers 
& Benson 1983) and low luminosity IRAS sources (Beichman et al. 1986). The magnetic 
field strengths (B) in the dense cores are estimated to be ~ 60 — 300/iG, assuming that 
the field within the cloud scales as B = Bo(n/no) K (Mouschovias 1978), where | < k < |, 
no = 1 cm -3 and Bo = 3/iG. The range given for k is governed by observational constraints 
and theoretical calculations for the lower and upper limits, respectively, the latter more 
closely defined than the former (see discussion in Mouschovias 1978). 

Dust grains in the outer layers of the cloud (Ay < 3) appear to have optical properties 
similar to those in the diffuse ISM (Vrba & Rydgren 1984; Kenyon et al. 1994). In lines of 
sight that sample regions of higher density, infrared spectroscopy of the 3.0/xm water-ice 
feature demonstrates the presence of icy mantles on the grains (Whittet et al. 1988). The 
transition from bare to mantled grains appears to occur at a rather well-defined value of 
Ay = 3.3 ± 0.1 (the threshold extinction). For lines of sight with extinctions higher than 
this threshold value, the optical depth of the 3.0/xm feature (and hence the column density 
of water molecules on the grains) increases linearly with Ay. The feature of solid CO at 
4.67/xm demonstrates a similar behavior, with a higher threshold extinction of Ay ~ 6 mag 
(Whittet et al. 1989, Chiar et al. 1995). Mantle growth represents a total change in grain 
surface properties, upon which some alignment mechanisms critically depend, and it is 
therefore of interest to compare p/A at extinctions above and below these thresholds. 

In the following sections, we investigate the correlation between p/A and A in the near 
infrared (A' band; Ao = 2.2/xm) for field stars background to the Taurus dark cloud. We 
make use of previous measurements of polarization together with a refinement of values 
for extinction from the literature. In § 2 we present our data compilation and results. 
Implications are discussed in § 3, and we summarize our conclusions in § 4. 
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2 THE DATA 



As compilations of visual polarization are heavily biased to low-extinction objects, 
observations in the K passband offer a better opportunity to construct a homogeneous set 
of data covering stars with a wide range of extinctions. Values of pk listed in Table 1 
are derived from Moneti et al. (1984), Tamura et al. (1987), Goodman et al. (1992) and 
Whittet et al. (1992, 1995). Data were rejected in cases where the fractional error (o' Pk /pk) 
is found to be greater than 20%. 

Measurements of pk and explicitly given in three references (Whittet et al. 

1995, Tamura et al. 1987 and Goodman et al. 1992). Polarimetric data from Moneti et al. 
(1984) are listed in the V passband and converted using an average ratio oipxjpv = 0.20, 
which was calculated using Taurus objects from Whittet et al. (1992). Note that the 
Pk values of Whittet et al. (1992) are measured in a non-standard passband centered at 
2.04/im. In these cases, our px values are determined using a power law extrapolation of 
the listed values of p at 1.64 and 2.04/xm (Whittet et al. 1992, Table 3), except in four 
cases where the value at 2.04/im was unknown (see Table 1). The Serkowski formula 

p/Pmax = exp[-A'ln 2 (A max /A)] (1) 

(Serkowski 1973; Coyne et al. 1974; Serkowski et al. 1975) was used to derive a value of 
Pk for these stars, with the parameters listed in Table 4 of Whittet et al. (1992). 

Corresponding extinction data (Ak) are also listed in Table 1. In the cases where the 
spectral type is known, visual and infrared photometry from the literature (see Table 1 for 
references) are used to derive infrared extinctions (Ak) assuming a standard extinction 
law: A K = O.lAy, where A v = I.IE V -k (Whittet k van Breda 1978). 

In the cases where the spectral types are not known (as in Tamura et al. 1987 and 
Goodman et al. 1992), a spectral type of K5III and an intrinsic color of (J — K)o = 0.95±0.4 
are assumed. This spectral type corresponds to the mean of the Elias (1978) set of objects, 
where the error in intrinsic color reflects the range in (J — K)o from types KOIII through 
M8III (the span of the Elias 1978 set). Intrinsic colors are taken from Bessell & Brett 
(1988). 

For stars with the highest extinctions, only photometry in the near-infrared passbands 
(J, H } K) could be found in the literature. Hence, we use an empirical formula derived 
by Whittet (1992) for extinction based upon infrared color excess, Ej-k- 

A k = 0.1 r Ej-k , (2) 
where r may be determined from the total to selective extinction ratio (Rv) by the formula 

2.332 



0.778 - l.lQAR' 1 



(3) 



A total of 14 Taurus sources whose photometric magnitudes in the B } V and K 
passbands are known (Whittet et al. 1995) give a mean value of Ry ~ 3.1, leading to 
r = 5.8 from equation (3). Thus, we use this value in equation (2) in our calculations of 
A K - 
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In the case of four stars from Whittet et al. (1995), we combine known values of 
A max and Eb-v with the relation Ry = 5.6 A max (Whittet & van Breda 1978) to derive 
Ak = 0.56 A max Eb-v (see Table 1). 

In Figure 1, we plot the resulting values of pk /Ak as a function of Ak for the stars 
listed in Table 1. 

3 DISCUSSION 

Figure 1 displays a monotonically decreasing relationship between alignment efficiency 
(p/A) and extinction (A), indicating that the linear polarization of radiation due to aligned 
dust grains in the Taurus dark cloud is produced more efficiently in regions of low extinction 
than in the more obscured lines of sight. Several physical processes may contribute to this 
trend in p/A. Similar decreases in polarization efficiency in the outer regions of dark clouds 
have been previously observed in p Ophiuchi (Vrba et al. 1976,1993), R Coronae Australis 
(Vrba et al. 1981), and Chamaeleon I (McGregor et al. 1994; Whittet et al. 1994). It may 
be worth noting that similar declines with extinction have been observed for the diffuse 
band (DIB) production efficiencies (W\/ Eb-v) hi p Oph by Snow & Cowen (1974) and 
in Taurus by Adamson et al. (1992), although no direct relationship between DIB's and 
aligned interstellar dust grains has ever been established. 

If one considers only the effects of density (n) on classical DG alignment, it is possible 
to derive the density dependence of the cloud's magnetic field strength, B (following the 
method of Vrba et al. 1981). For grains of equal radius, 

P/A pg = 2mFco^u { ( aa/ab) - 1 , (4) 

(Davis & Greenstein 1951), where P is the polarization in magnitudes, A pg is the pho- 
tographic extinction (A pg Ab), F is the distribution parameter of grain orientations, 
(■| — v) is the angle between the direction of the magnetic field (B) and the line of sight, 
and a a and cr& are the extinction coefficients parallel and perpendicular to the long axes of 
the grains, respectively. Assuming that B is perpendicular to the line of sight, then cos 2 v 
= 1. If we take o^j ' o a = b/a (Davis & Greenstein 1951, Vrba et al. 1981) and a typical in- 
terstellar grain axial ratio of b/a = 0.2 (Aanestad & Purcell 1973), equation (4) reduces to 
P/A pg = 1.2 F. Conversion of A pg to A K and P to p (in %) gives p/A K = 730 F% mag" 1 . 
From Jones & Spitzer (1967), 

^ ? = ^ 1 ( J ^) 1/2 (7-l)(l-% L ), (5) 
(baton mkl 1 

where the imaginary part of the grain's magnetic susceptibility (x") and its rotational 
velocity (u) are related as = 2.6 X lO -12 !^" 1 for reasonable grain materials and values 

of to (Purcell 1969). T gr is the grain temperature, T the surrounding gas temperature, B 
the strength of the magnetic field, a the grain size, and 7 the ratio of the moments of 
inertia about the short and long axes, respectively, 7 = ^[(a/b) 2 + 1] (Davis & Greenstein 
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1951). Taking T = 50K, T gr = 10K as typical for dense clouds outside of dense cores and 
combining equations (4) and (5), we find 

5.7 x 10 8 B 2 

P/Ak = % mag . (6) 

an 

Assuming a density dependence for the magnetic field strength of B = Bo(n/no) K (Mousch- 
ovias 1978) and using the mean interstellar values of Bo = 3 /iG and no = 1cm -3 , we derive 

p /A K = 5 - lxl0 ~ 3 ra 2«-i % -i. (7) 

a 

Since n = N / L (where N is the column density of the cloud and L is the length of the 
line of sight through it) and it is known from observations that A oc iV(H), the power law 
dependence on n in equation (7) also holds for A. Thus, p/Ak oc A 2 ^ -1 ■ 
A least-squares power law fit to all data points in Table 1 yields 

Pk/Ak = 1-38 A" ' 53 % mag" 1 . (8) 

Relating this to equation (7), we derive a scaling factor for B of n = 0.24 for the Taurus 
dark cloud. 

This model considers the scaling of B with n, ignoring the effect of grain growth. The 
relative enhancements of the ratio of total to selective extinction (Rv) and the wavelength 
of maximum polarization (A max ) toward dense regions of dark clouds suggest that grain 
growth by mantle condensation or grain-grain coagulation occurs there, producing grains 
which are more spherical and thus reducing the efficiency with which they may be aligned 
(Carrasco et al. 1973; Whittet & van Breda 1975; Whittet & Blades 1980). Models of 
polarization due to alignment by the DG mechanism show that in fact grains which are 
less elongated do have lower values of p/A (Purcell & Spitzer 1971; Aannestad & Purcell 
1973). On the basis of the variation in A max with Ay in the p Oph cloud, Vrba et al. (1993) 

1/12 

deduce that the effect of grain growth may be described as a = 1.18 a{ Ay , where a{ is 
the initial value of a. This introduces an additional factor of n -1 / 12 to equation (7), and 
hence, taking this into account, k = 0.28. This value is somewhat below the lower limit of 
| given by Mouschovias (1978), but within a reasonable range considering the assignment 
of this limit, and also consistent with the results of Tamura et al. (1987) and Goodman et 
al. (1992). 

Stronger declines in p/A with extinction may occur due to several different processes 
within dense interstellar clouds. Grain surfaces are active sites for the formation of H2. A 
hydrogen atom colliding with a dust grain has a high probability of sticking, followed by 
migration over the surface until finally combining with another hydrogen atom (e.g. Tielens 
& Allamandola 1987). The energy of this formation process may be transferred in part to 
kinetic energy, ejecting the molecule from the surface and donating angular momentum to 
the grain. Molecular formation is likely to occur at preferential sites on the grain surface 
due to defects or impurities in the surface structure (Hollenbach & Salpeter 1971), and if 
the number of active sites are small, then their distribution over the surface will determine 
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the spin properties of the grain (Purcell 1975, 1979). A series of recombination events 
occurring at these preferential sites could then lead to grain angular velocities that are a 
factor of ~ 10 5 higher than those predicted by random elastic gas-grain collisions. This 
suprathermal spin effect would lead to very efficient alignment of interstellar grains, but 
the physical processes required may not be available in all grain environments. 

Undoubtedly, the largest influence on the suprathermal spin mechanism is that the 
abundance of atomic hydrogen available with which to form H2 falls rapidly with increasing 
cloud density (Savage et al. 1977). Moreover, declines in p/A could be expected solely on 
the basis of H2O ice mantle growth since hydrogen adsorption sites on an H2O surface 
are weak physisorption sites, as opposed to the stronger sites obtainable on a bare silicate 
or carbonaceous surface. Hence, the time a hydrogen atom spends on an H2O surface 
(and thus the probability that it will combine with a second such adsorbed H atom) will 
be smaller than that on a bare, unmantled grain (Tielens & Allamandola 1987, Williams 
1993). 

Mathis (1986) suggest that a small inclusion of superparamagnetic materials (e.g. 
metallic Fe, Fe3 0"4, etc) in interstellar grains could increase the imaginary part of their 
magnetic susceptibility (x") by a factor of 10 6 over ordinary materials. While super- 
paramagnetic (SPM) alignment could efficiently align grains even in the observed galactic 
magnetic field, its environment dependence is not different from that of DG alignment. 

4 CONCLUSIONS 

In combining several data sets for the Taurus dark cloud for a large range in extinction 
(Ak < 2.5 mag, Ay < 25 mag), we have shown that the efficiency with which background 
starlight is polarized by dust (p/A) decreases strongly with the amount of obscuring ma- 
terial in the line of sight. As previously noted, the form of Fig 1 is a monotonic decrease 
in p/A with A. We have shown that this behavior is consistent with the expected depen- 
dence of B on n for DG-type alignment. Changes in grain surface properties associated 
with the appearance of H2O mantles (Ak > 0.3) and CO mantles (Ak > 0.6) do not 
have an obvious effect on the alignment efficiency within the scatter of our data, in spite 
of the fact that grains with icy surfaces should possess different suprathermal spin prop- 
erties than bare grain cores. The presence (or absence) of mantles may have little effect 
on the ability of a grain to align in the ambient magnetic field. This is consistent with 
the observed presence of polarization enhancement in the 3/xm ice feature (Hough et al. 
1988, 1989) which demonstrates that grains mantled with water-ice do indeed align. Small 
but non-zero alignment efficiencies deep within dark clouds are a natural consequence of 
alignment by the DG mechanism. 
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Table 1. Spectral type, Aj{ } Pk, °> k ■> an d Pk/Ak for the set of objects studied 
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Table 1. continued 



Star I.D. Sp Type a Ak Pk Cp K Pk/Ak Ref 

mag % % % mag -1 

RA 4 21 27.8 — 0.14 1.29 0.09 9.21 16 

RA 4 21 29.3 — 0.14 0.65 0.08 4.64 16 

RA 4 22 11.6 — 0.55 1.01 0.20 1.84 16 

RA 4 22 34.4 — 0.43 1.47 0.26 3.42 16 

Tam 8 K5III 2.33 2.7 0.3 1.16 4,17 

Tam 12 — 0.66 1.5 0.2 2.27 4 

Tam 16 — 0.60 2.0 0.2 3.33 4 

Tam 17 — 0.60 1.2 0.1 2.00 4 



Notes - (a) If no spectral type is listed, type K5III has been assumed. In these cases, 
values of Ak contain an error of about 0.2 mag due to uncertainty in assigned 
intrinsic color, (b) pk estimated using the values of p max , Vax and K from Whittet 
et al. (1992). Fractional error (& Pk /pk) taken equal to that of Cp max /pmax- 
(c) A K = 0.56A max E B -v- 

References- (1) Elias 1978; (2) Whittet et al. 1995; (3) Moneti et al. 1984; 
(4) Tamura et al. 1987; (5) Whittet et al. 1992; (6) Straizys & Meistas 1980; 
(7) Vrba k Rydgren 1985; (8) Kenyon et al. 1994; (9) Oja 1987; 
(10) Crutcher 1985; (11) Straizys et al. 1982; (12) Simbad database; 
(13) Ungerer et al. 1985; (14) Straizys 1985; (15) Slutskii et al. 1980; 
(16) Goodman et al. 1992; (17) Whittet et al. 1988. 
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Figure 1. Values of pk/A-k plotted against Ax (from Table 1). Filled circles indicate 
the stars with known spectral types. Stars with uncertain spectral types are plotted 
using empty circles. The dotted line represents a least-squares fit to all data points: 
p K /A K = 1.38V 53 - 
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